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Abstract 

The relativistic proton-neutron quasiparticle random phase approximation (PN-RQRPA) is ap- 
phed in the calculation of total muon capture rates on a large set of nuclei from ^^C to ^^^Pu, for 
which experimental values are available. The microscopic theoretical framework is based on the 
Relativistic Hartree-Bogoliubov (RHB) model for the nuclear ground state, and transitions to ex- 
cited states are calculated using the PN-RQRPA. The calculation is fully consistent, i.e., the same 
interactions are used both in the RHB equations that determine the quasiparticle basis, and in the 
matrix equations of the PN-RQRPA. The calculated capture rates are sensitive to the in-medium 
quenching of the axial-vector coupling constant. By reducing this constant from its free-nucleon 
value gA = 1.262 by 10% for all multipole transitions, the calculation reproduces the experimental 
muon capture rates to better than 10% accuracy. 

PACS numbers: 21.60.Jz, 23.40.Bw, 24.30.Cz, 25.30.Mr 
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I. INTRODUCTION 



Semi-leptonic weak interaction processes in nuclei are very sensitive to detailed properties 
of nuclear ground states and excitations. In astrophysical applications, in particular, weak 
interaction rates (/?-decay half lives, neutrino-nucleus cross sections, electron capture rates) 
must be calculated for hundreds of isotopes. Many of those are located far from the valley 
of /3-stability, and thus not easily accessible in experiments. For a consistent description, 
reliable predictions and extrapolations of these processes it is, therefore, essential to employ 
a consistent theoretical framework based on microscopic nuclear structure models. 

At present the framework of nuclear energy density functionals (NEDF) provides the 
most complete description of ground-state properties and collective excitations over the 
whole nuclide chart. At the level of practical applications the NEDF framework is realized in 
terms of self-consistent mean-field (SCMF) models. With a small set of universal parameters 
adjusted to data, the SCMF approach has achieved a high level of accuracy in the description 
of structure properties over the whole chart of nuclides, from relatively light systems to 
superheavy nuclei, and from the valley of /5-stability to the particle drip- lines {l|, 3] ■ 

In a series of recent studies we have used a fully consistent microscopic approach based on 
relativistic energy density functionals to analyze /?-decay half-lives of neutron-rich nuclei 3|, 
^, and to model inclusive charged-current neutrino- nucleus reactions [5]. In this framework 
nuclear ground states are described using the relativistic Hartree-Bogoliubov (RHB) model 
l2i, and transitions to excited nuclear states are calculated in the relativistic quasiparticle 

f] Q 

random-phase approximation (RQRPA) [6|, |7[. There are important advantages in using 
functionals with manifest covariance, the most obvious being the natural inclusion of the 
nucleon spin degree of freedom. The resulting nuclear spin-orbit potential has the correct 
empirical strength and isospin dependence. This is, of course, especially important in the 
description of excitations in the spin-isospin channel, e.g. semi-leptonic weak interaction 
processes. In addition, by employing a single universal effective interaction in modeling 
both ground-state properties and multipole excitations in various mass regions of the chart 
of nuclides, the calculation of weak-interaction rates is essentially parameter free, and can 
be extended to regions of nuclei far from stability, including those on the r-process path. 

To successfully extend a particular microscopic approach to regions of unknown nuclei far 
from stability, it is necessary to perform extensive tests and compare results with available 
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data. Reliable prediction of weak interaction rates, in particular, require a fully consistent 
description of the structure of ground states and multipole excitations. For instance, calcu- 
lated /3-decay half-lives are very sensitive to low-energy Gamow- Teller transitions, but can 
only test excitations of lowest multipoles. Higher multipoles are excited in neutrino-nucleus 
reactions in the low-energy range below 100 MeV, and these reactions could play an im- 
portant role in many astrophysical processes, including stellar nucleosynthesis. There are, 
however, only few data on neutrino-nucleus reactions, and these are limited to relatively 
light nuclei. Much more data are available for total muon capture rates. Muon capture on 
stable nuclei has been studied in details since many years, both experimentally and theo- 



retically 



0, y, Q, 



In this process the momentum transfer is of the order of the muon 



mass and, therefore, the calculation of total muon capture rates presents an excellent test 
of models that are also used in studies of low-energy neutrino-nucleus reactions. 

In this work we test the fully consistent RHB plus proton-neutron RQRPA model in 
the calculation of total muon capture rates on a large set of nuclei from ^^C to ^'^'^Pu, for 



which experimental values are available y^. Previous calculation of muon capture rates 
on selected nuclei using the RPA approach include the consistent Hartree-Fock (HF) RPA 

n n 

model [13|, [IJ], in which the HF mean field and the particle- hole interaction result from 
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18| in which both the 



the same Skyrme effective force, and a series of studies 

continuum and standard RPA were used, and the effect of quenching of axial- vector coupling 

was analyzed. The present analysis parallels the recent study by Zinner, Langanke and Vogel 

jisl ] , where the nonrelativistic RPA was used to systematically calculate muon capture rates 

for nuclei with 6 < Z < 94. There are, however, significant differences between the two 
- - ^ 

approaches. The model employed in Refs. [15|, 



16 



13, 



18j uses a phenomenological Woods- 



Saxon potential to generate the basis of single-nucleon states. The strength of the potential 
is adjusted to experimental proton and neutron separation energies in individual nuclei. In a 
second step the RPA with a phenomenological Landau-Migdal residual interaction is used to 
calculate nuclear excitations. The present approach, as already emphasized above, is fully 
consistent: both the basis of single-nucleon states and multipole excitations of nuclei are 
calculated from the same energy density functional or nuclear effective interaction. Results 
will be compared with data and discussed in relation to those reported in Ref. ISj. In 
particular, we will consider the important issue of quenching of the axial-vector strength. 
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II. THEORETICAL FRAMEWORK 



The capture of a negative muon from the atomic Is orbit on a nucleus {Z, N) 

Ij- + {Z,N) ^u^ + {Z -i,N+iy , (1) 

presents a simple semi-leptonic reaction that proceeds via the charged current of the weak 
interaction. Detailed expressions for the reaction rates and the transition matrix elements 



can be found in Refs. 



11 



19 



20|. The capture rate reads 



lepton spins ^^/j "h j- ^j^. 



w V 



(2) 



where fl denotes the quantization volume and u is the muon neutrino energy. The Hamil- 
tonian Hw of the weak interaction is expressed in the standard current-current form, i.e. in 
terms of the nucleon Jx{x) and lepton i\{x) currents 



G f 

w = J dxJx{x)i^{x) 



and the transition matrix elements read 



G 



{f\Hw\i) = -^/a 



1/v^ 



qx 



(3) 



(4) 



(f)is{x) is the muon Is wave function, the four-momentum transfer is g = {qo,q), and the 
multipole expansion of the leptonic matrix element Le"^'^'^ determines the operator structure 



for the nuclear transition matrix elements 
rate is given by 

2GV2_ 1 ' 

. j=o 



11 
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20|]. The expression for the muon capture 



UJfi 



;i + u/Mt) 2Ji + 



'Is 



Mj - C 



J. 



j=i 



J^ 

(5) 



where G is the weak coupling constant, the phase-space factor (1 + v/Mt)~^ accounts for 
the nuclear recoil, and Mt is the mass of the target nucleus. The nuclear transition ma- 
trix elements between the initial state \ Ji) and final state |J/), correspond to the charge 
yWj, longitudinal £j, transverse electric Tf^^ and transverse magnetic Tj'^^'^ multipole 
operators: 



the Coulomb operator 



MjM{x) = F^Mf{x)~i 



K 



rriM 



FAQfix) + -{Fa - m,Fp)i:f\x) 



(6) 
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the longitudinal operator 



(7) 



the transverse electric operator 



JM 



K 



TUn 



(8) 



and the transverse magnetic operator 



K 



+ FAJ:fix 



(9) 



where all the form factors are functions of g^, and n = \q\. These multipole operators contain 
seven basic operators expressed in terms of spherical Bessel functions, spherical harmonics, 
and vector spherical harmonics 19 L By assuming conserved vector current (CVC), the 
standard set of form factors reads 2l| : 



/i^(g2) = 4.706 

Fp{q 



1 + 



840 MeV 



1 -2 



-1.262 



1+1 ^ V 
' 840 MeV J 



1032 MeV 



2^ _ 2mAfFA(g^) 

g2 _|_ ^2 



(10) 

(11) 

(12) 
(13) 



The muon capture rates are evaluated using Eq. ([5]) , with the transition matrix elements 
between the initial and final states determined in a fully microscopic theoretical framework 
based on the Relativistic Hartree-Bogoliubov (RHB) model for the nuclear ground state, 
and excited states are calculated using the relativistic quasiparticle random phase approx- 
imation (RQRPA). The RQRPA has been formulated in the canonical single-nucleon basis 



of the relativistic Hartree-Bogoliubov (RHB) model in Ref. {g], and extended to the descrip- 
tion of charge-exchange excitations (proton-neutron RQRPA) in Ref. In addition to 
configurations built from two-quasiparticle states of positive energy, the relativistic QRPA 
configuration space must also include pair-configurations formed from the fully or partially 



occupied states of positive energy and empty negative-energy states from the Dirac sea. 
The RHB+RQRPA model is fully consistent: in the particle-hole {ph) channel effective La- 
grangians with density-dependent meson-nucleon couplings are employed, and pairing (pp) 
correlations are described by the pairing part of the finite range Gogny interaction. Both 
in the ph and pp channels, the same interactions are used in the RHB equations that deter- 
mine the canonical quasiparticle basis, and in the matrix equations of the RQRPA. In this 
work we use one of the most accurate meson-exchange density-dependent relativistic mean- 
field effective interactions - DD-ME2 {2^ in the ph channel, and the finite range Gogny 
interaction DIS jQ] in the pp channel. 

The spin-isospin-dependent interaction terms are generated by the vr- and p-meson ex- 
change. Although the direct one-pion contribution to the nuclear ground state vanishes at 
the mean-field level because of parity conservation, the pion must be included in the calcu- 
lation of spin-isospin excitations. The particle-hole residual interaction of the PN-RQRPA 
is derived from the Lagrangian density: 

4^*, = -g^ijjf^p^ri; - A^^^^M^^^f^ (14) 

where vectors in isospin space are denoted by arrows. For the density-dependent coupling 
strength of the p-meson to the nucleon we choose the value that is used in the DD-ME2 

n 

effective interaction 22!], and the standard value for the pseudovector pion-nucleon coupling 
is /^/47r = 0.08, and = 138 MeV. The derivative type of the pion-nucleon coupling 
necessitates the inclusion of the zero-range Landau-Migdal term, which accounts for the 
contact part of the nucleon-nucleon interaction 
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V,. = g'\^] rir25]i-S2(5(r-i-r2), (15) 



with the parameter g' adjusted in such a way that the PN-RQRPA reproduces experimental 
values of Gamow- Teller resonance (GTR) excitation energies [3]. The precise value depends 
on the choice of the nuclear symmetry energy at saturation, and for the DD-ME2 effective 
interaction (7'=0.52 has been adjusted to the position of the GTR in ^''^Pb. This value is 
kept constant for all nuclides calculated in this work. 

In the evaluation of muon capture rates (Eq. ([5])), for each transition operator Oj the 
matrix elements between the ground state of the even-even (A^, Z) target nucleus and the 
final state are expressed in terms of single-particle matrix elements between quasiparticle 



canonical states, the corresponding occupation probabilities and RQRPA amplitudes: 



{Jf\\dj\\Ji) = Y.iP\\0'^\\^) [X^pnUpVn - Ypn^pUn 



(16) 



pn 



Transitions between the |0+) ground state of a spherical even-even target nucleus and excited 
states in the corresponding odd-odd nucleus are considered. The total muon capture rate is 
calculated from the expression: 



2gH y - 

1 f^n I 



j% (1 + ^//Mt) 



0^ 



Jf 



>ls 



f rj-el rj-mag 



j% (1 + ^s/Mt) 

with the neutrino energy determined by the energy conservation relation 



(17) 



'mf,-eb + Ei = Ef + Uf 



(18) 



where is the binding energy of the muonic atom. 

For each nucleus the muon wave function and binding energy are calculated as solu- 
tions of the Dirac equation with the Coulomb potential determined by the self-consistent 
ground-state charge density. However, while the RHB single-nucleon equations are solved 
by expanding nucleon spinors and meson fields in terms of eigenf unctions of a spherically 
symmetric harmonic oscillator potential, the same method could not be used for the muon 
wave functions. The reason, of course, is that the muon wave functions extend far beyond 
the surface of the nucleus and, even using a large number of oscillator shells, solutions ex- 
pressed in terms of harmonic oscillator basis functions do not converge. The Dirac equation 
for the muon is therefore solved in coordinate space using the method of finite elements with 
B-spline shape functions [24, l25|. As an illustration, in Fig. [T]we plot the square of the Is 
muon wave functions for ^^O, ^°Ca, ^^°Sn and ^''^Pb. The solutions that correspond to self- 
consistent ground-state charge densities are compared with eigenf unctions of the Coulomb 
potential for the corresponding point-charge Z. For light nuclei the radial dependence of the 
Is muon wave function is not very different from that of the point-charge Coulomb potential. 
With the increase of Z the muon is pulled into the nuclear Coulomb potential, and thus the 
magnitude of the Is density inside the nucleus is reduced with respect to the point-charge 
value. To test our calculation of muon orbitals in the nuclear Coulomb potential, in Tables [T] 
andlUthe muon transition energies in Sn isotopes and in ^°^Pb, respectively, are compared 



with available data 26|, |27|] . The calculated transition energies are in good agreement with 
experimental values. 

The effect of the finite distribution of ground-state charge densities on the calculated 
muon capture rates is illustrated in Fig. [21 For a large set of nuclei from ^^C to ^'^^Pu, we 
plot the ratio between calculated and experimental muon capture rates. This ratio is < 1.5 
for all nuclei when the muon Is wave functions are determined by self-consistent ground- 
state charge densities, whereas for point-charge Coulomb potentials one notes a distinct 
increase with Z, and cUcaic./i^exp. > 4 for the heaviest systems. 

III. RESULTS FOR MUON CAPTURE RATES 

The muon capture rates shown in Fig. [2] are calculated with the standard set of free 
nucleon weak form factors Eqs. fllOp - (1131) 2l| . i.e. the calculation does not include any in- 



medium quenching of the corresponding strength functions. Even with muon wave functions 
determined self-consistently by finite-charge densities, the resulting capture rates are larger 
than the corresponding experimental values by a factor ^ 1.2 — 1.4. This is in contrast 



to the results of Ref. 18|], where the experimental values have been reproduced to better 
than 15% accuracy, using the free-nucleon weak form factors and residual interactions with 
a mild A dependency. In fact, it was shown that the calculated rates for the same residual 
interactions would be significantly below the data if the in-medium quenching of the axial- 
vector coupling constant is employed to other than the true Gamow- Teller (GT) amplitudes. 



Consequently, the calculations reported in Ref. 18| were performed with quenching only the 
GT part of the transition strength by a common factor (0.8)^ = 0.64. It was concluded, 
however, that there is actually no need to apply any quenching to operators that contribute 
to the muon capture process, especially those involving single-nucleon transitions between 
major oscillator shells. 

As already emphasized in the Introduction, although both calculations are based on the 
RPA framework, there are important differences between the model of Ref. 
RHB+RQRPA approach employed in the present study. The main difference is probably 
the fact that the present calculation is fully consistent: for all nuclei both the basis of single- 
nucleon states and the multipole response are calculated using the same effective interaction, 
whereas in Ref. jl8] the phenomenological Woods-Saxon potential was adjusted to individual 



18| . and the 
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nuclei and the strength of the residual Landau-Migdal force had a mild A-dependence. 

In Fig. [3] we compare the ratios of the theoretical and experimental total muon capture 
rates for two sets of weak form factors. First, the rates calculated with the free nucleon 
weak form factors Eqs. flTU]) - f[T^ 21| (circles), and already shown in Fig. [21 The lower 
rates, denoted by diamonds, are calculated by applying the same quenching qa = 1.262 — > 
qa = 1.135 to all axial operators, i.e. qa is reduced by 10% in all multipole channels. 
In the latter case the level of agreement is very good, with the mean deviation between 
theoretical and experimental values of only 6%. The factor 0.9 with which the free-nucleon 
qa is multiplied is chosen in such a way to minimize the deviation from experimental values 
for spherical, closed-shell medium-heavy and heavy nuclei. On the average the results are 



slightly better than those obtained in Ref. [18 
calculation of Zinner, Langanke and Vogel 



cf. Fig. 2 of [18|). Note, however, that in the 



18l | only the true Gamow- Teller Qhuj transition 



strength was quenched, rather than the total strength in the l"*" channel. In the present 
study considerably better results are obtained when the quenched value of the axial-vector 
coupling constants is used for all multipole operators. The reason to consider quenching the 
strength in all multipole channels, rather than just for the GT is, of course, that the axial 
form factor appears in all four operators Eqs. ([6]) - ([9]) that induce transitions between the 
initial and final states, irrespective of their multipolarity. Even more importantly, only a 
relatively small contribution to the total capture rates actually comes from the GT channel 
1+. This is illustrated in Fig. HI where we display the relative contributions of different 
multipole transitions to the RHB plus RQRPA muon capture rates in ^^O, '^°Ca, ^^°Sn and 
^°^Pb. For the two lighter N = Z nuclei the dominant multipole transitions are A'^ = 1~ and 
A'^ = 2^ (spin-dipole). For the two heavier nuclei there are also significant contributions of 
the A'^ = l"*" and A'^ = 2+, especially for ^osp]^ ^^^^ f^j, other heavy nuclei. Note that in heavy 
nuclei the A'^ = 1^ multipole represents 2huj transitions, rather than the Ohu Gamow- Teller 
transitions. 

Returning to Fig. [3l we notice that with a 10% quenching of the free-nucleon axial- 
vector coupling constant qa, for medium-heavy and heavy nuclei the calculated capture 
rates are still slightly larger than the corresponding experimental values, with the ratio 
locale. /i^exp. typically around 1.1, whereas for several lighter nuclei considered here this ratio 
is actually less than 1 (cf. also Table HTB . Overall the best results, with ci^caic./i^exp. ~ 1, are 
obtained near closed shells. The characteristic arches between closed shells can probably 
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be attributed to deformation effects, not taken into account in our RHB+RQRPA model. 
In addition to the DD-ME2 interaction, we liave also carried out a full calculation of total 
capture rates from ^^C to ^'''^Pu, using the density- and momentum-dependent relativistic 
effective interaction DSC*. In the study of /5-decay half-lives of Ref. ^, this interaction was 
constructed with the aim to enhance the effective (Landau) nucleon mass, and thus improve 
the RQRPA description of /5-decay rates. When DSC* is used to calculate muon capture 
rates, some improvement is obtained only locally, for certain regions of Z, whereas in other 
regions {Z ^ 50 and Z > 82) the results are not as good as those obtained with DD-ME2. 
The overall quality of the agreement between theoretical and experimental capture rates is 
slightly better with DD-ME2. 

The calculated total muon capture rates for natural elements and individual isotopes 
are also collected in Table UTTl and compared with available data 12|. In particular, the 
calculation nicely reproduces the empirical isotopic dependence of the capture rates Q], i.e. 
for a given proton number Z the rates decrease with increasing neutron number, because 
of the gradual blocking of available neutron levels. The isotopic trend is also illustrated in 
Fig. [5l where we plot the experimental and theoretical total muon capture rates on Ca, Cr 
and Ni nuclei. The latter correspond to the quenching qa = 1.262 —>■ Qa = 1-1S5 for all 
multipole operators. 

In conclusion, we have tested the RHB plus proton-neutron RQRPA model in the calcula- 
tion of total muon capture rates on a large set of nuclei from ^^C to ^^"^Pu. The calculation is 
fully consistent, the same universal effective interactions are used both in the RHB equations 
that determine the quasiparticle basis, and in the matrix equations of the RQRPA. The cal- 
culated capture rates are sensitive to the in-medium quenching of the axial-vector coupling 
constant. By reducing this constant from its free- nucleon value qa = 1.262 to the effective 
value Qa = 1.1S5 for all multipole transitions, i.e. with a quenching of approximately 10%, 
the experimental muon capture rates are reproduced with an accuracy better than 10%. 
This result can be compared to recent RPA-based calculations Q, 17, 18|, that reproduce 
the experimental values to better than 15%, using phenomenological potentials adjusted to 
individual nuclei and A-dependent residual interactions, but without applying any quench- 
ing to the operators responsible for the fi~ capture process. The test has demonstrated that 
the RHB plus QRPA model provides a consistent and accurate description of semi-leptonic 
weak interaction processes at finite momentum transfer in medium-heavy and heavy nuclei 
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over a large Z-range. The fully consistent microscopic approach, based on modern relativistic 
nuclear energy density functionals, can be extended to other types of weak interaction pro- 
cesses (electron capture, neutrino- nucleus charge-exchange and neutral-current reactions), 
and to regions of short-lived nuclei far from stability. 
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FIG. 1: The square of the Is muon wave function in the Coulomb potentials of self-consistent 
ground-state charge densities of ^^O, ^'^Ca, ^^'^Sn and 208p|^ (solid curves), compared to eigenfunc- 
tions of the Coulomb potential for the corresponding point charge Z (dashed curves). The figures 
also include the calculated charge densities of the four nuclei, scaled by arbitrary factors. 



13 




FIG. 2: Ratio of the calculated and experimental total muon capture rates, as function of the 
proton number Z. The theoretical values are calculated with muon Is wave functions determined 
by self-consistent ground-state charge densities (filled circle symbols), and by the corresponding 
point-charge Coulomb potentials (squares). 
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FIG. 3: (Color online) Ratio of the calculated and experimental total muon capture rates, as 
function of the proton number Z. Circles correspond to rates calculated with the free-nucleon 
weak form factors Eqs. (jlOp - (|13p 



21[, and diamonds denote values obtained by quenching the 
free-nucleon axial-vector coupling constant = 1-262 to gA = 1.135 for all operators, i.e. in all 
multipole channels. 
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FIG. 4: (Color online) Relative contributions of different multipole transitions to the RHB plus 
RQRPA total muon capture rates in ^"^O, ^°Ca, ^^OSn and ^ospb. 
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FIG. 5: (Color online) Total muon capture rates on Ca, Cr and Ni isotopes. Experimental rates 
(filled symbols) are compared to theoretical values (empty symbols), calculated using the fully 
consistent RHB plus RQRPA framework with the DD-ME2 universal effective interaction, and 
with the quenching of the axial-vector coupling constant = 1-262 ^ qa = 1-135 for all multipole 
operators- 
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TABLE I: Calculated muon transition energies in tin isotopes (in units of keV), compared with 



available data 
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i20Sn 


3408 


3415 


3454 


3460 


i22Sn 




3409 




3454 


i24Sn 


3400 


3404 


3445 


3450 



TABLE II: Calculated muon transition energies in ^°^Pb (in units of keV), in comparison with 
experimental values 



208 


Pb 


exp. 


calc. 


lP3/2 - 


- l^i/a 


5963 


5956 


lPl/2 - 


- l'5l/2 


5778 


5773 


1^3/2 - 


- lPl/2 


2642 


2633 


1^5/2 - 


- lP3/2 


2501 


2493 


lc^3/2 - 


- lP3/2 


2458 


2450 
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TABLE III: Experimental and calculated muon capture rates for natural elements and individual 
isotopes. The theoretical rates are calculated using the fully consistent RHB plus RQRPA frame- 
work with the DD-ME2 universal effective interaction, and with the quenching of the axial-vector 
coupling constant qa = 1.262 qa = 1.135 for all multipole operators. Values for naturally 
occuring elements (element symbol with no superscript) are weighted averages of capture rates on 
individual isotopes, using their natural abundances. Experimental values are from Ref. |l^, unless 
otherwise stated. All rates are in units of lO^s-^ 



Nucleus 


Exp. 


Calc. 


Nucleus 


Exp. 


Calc. 


Nucleus 


Exp. 


Calc. 


Nucleus 


Exp. 


Calc. 


12c 


0.039 





032 


78Se 




6.644 


i22Sn 




9.645 


164Dy 




13 


540 


160 


0.103 





065 


80Se 




5.796 


i24Sn 




8.837 


Dy 


12.29 


14 


194 


180 


0.088 





057 


82Se 




4.935 


Sn 


10.44 


10.923 


166gj, 




16 


129 


20Ne 


0.204 





237 


Se 


5.681 


5.950 


i26xe 




10.652 


168 




14 


949 


24Mg 


0.484 





506 


86Sr 




8.885 


i28xe 




9.830 


170Er 




13 


912 


28Si 


0.871 





789 


88Sr 




7.393 


isoxe 




9.068 


Er 


13.04 


15 


270 


32s 


1.352 


1 


485 


Sr 


7.020 


7.553 


Te 


9.270 


9.706 


ITSjjf 




16 


434 


40 


1.355 


1 


368 


90 Zr 




9.874 


i32xe 


9.4'' 


10.631 


180Hf 




15 


276 


40Ca 


2.557 


2 


340 


92 Zr 




9.694 


i36Xe 


8.6^ 


8.625 


Hf 


13.03 


15 


783 


^^Ca 


1.793 


1 


851 


94 Zr 




8.792 


i36Ba 




11.461 


182^ 




17 


259 


48Ca 


1.214'' 


1 
1 




Zr 


8.660 


9.619 


i38Ba 




10.127 


184w 








48xi 


2.590 


2 


544 


92Mo 




12.374 


Ba 


9.940 


10.259 


186w 




14 


807 


50Cr 


3.825 


4 


001 


94Mo 




12.001 


i40Ce 




11.888 


W 


12.36 


15 


971 


52Cr 


3.452 


3 


419 


96 Mo 




10.933 


i42Ce 




12.142 


198 jjg 




17 


369 


54Cr 


3.057 


3 


065 


98Mo 




9.804 


Ce 


11.60 


11.917 


200 jjg 




16 


227 


Cr 


3.472 


3 


483 


Mo 


9.614 


10.995 


142Nd 




14.043 


202Hg 




15 


205 




4.411 


4 


723 


104p^ 




13.182 






14.288 


204Hg 




13 


993 


58Ni 


6.110 


6 


556 


106pd 




11.912 


146Nd 




12.981 


Hg 


12.74 


15 


733 


60Ni 


5.560 


5 


610 


108pd 




10.795 


Nd 


12.50 


13.861 


206p]2 




15 


717 


62Ni 


4.720 


4 


701 






9.821 


i48Sm 




15.425 


208p]2 




13 


718 


Ni 


5.932 


6 


234 


Pd 


10.00 


11.391 


isosm 




14.132 


Pb 


13.45 


14 


348 


64Zn 




6 


862 


iiOCd 




12.960 


i52Sm 




13.451 


232rpj^ 


12.56 


13 


092 


66 Zn 




5 


809 


ii^Cd 




19 
11.800 


154c 
Sm 




12.563 


234u 


13.79 


14 


231 


68r7„ 




A 




114r^j 




-ir\ TAa 




1 o oo 


1 Q 


236 TT 


1 Q nnc 


1 Q 


Ann 



